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ABSTRACT 



We present deep IRAC images that highlight the scattered light emission 
around many of the youngest protostars, the so-called Class sources, in 
L1448. By comparison of the data with a Monte Carlo radiative transfer code 



(jWhitney et al.l l2003bl ) . we demonstrate for the first time that the observed in- 
frared light from these objects is consistent with scattered light from the central 
protostar. The scattered light escapes out the cavity, carved by molecular out- 
flows, in the circumstellar envelope. In particular, we observe prominent scattered 
light nebulae associated with the Class sources: L1448-mm, L1448 IRS 2, and 
3B, as well as a Class I source: IRS 3A. We use a grid of models with probable 
protostellar properties to generate model spectral energy distributions (SEDs) 
and images for bands sensitive to this scattered light: J, H, Ks, and Spitzer 
IRAC bands. By simultaneously fitting SEDs and images of the outflow cavities, 
we are able to model geometric parameters, i.e. inclination angle and opening 
angle, and loosely constrain physical parameters. The opening angle may be an 
important indicator of the evolutionary state of a source. We compare our results 
for Class sources to similar studies of Class I sources. There may be a transition 
phase from Class to Class I when a source has an opening angle between 20° to 
30°. It is important to note that while the best fit model image and SED do not 
fully describe the sources, the fits generally describe the circumstellar structure 
of Class sources in L1448. 

Subject headings: ISM: individual (L1448) — ISM: jets and outflows — stars: 
circumstellar matter — stars: formation — stars: pre-main sequence 
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Introduction 



The L1448 dark cloud is locate d on the western s ide of the Perseus molecular cloud com- 
plex at a distance of ~250 pc (e.g. Enoch et al.ll2006l ). The relatively isolated, young stellar 



population and their associated spectacular molecular outflows make L1448 an ideal star 



form ation laboratory. L1448 contains several embedded YSOs at the Class (lAndre et al. 



19931 ) and Class I stage, identified as such from infrared to millimeter observations. In this 



pape r 



1990 



, we study four sources i n particular: L1448-mm dCuriel et al. 



199ll. Il995h . L1448 IRS 2 flO'Linger et a 



20061), and L1448 IRS 3A and 3B flCuriel et al. 



1999 



Wolf- Chase et al. 



19901 ; lLooney et al. 



1990 



2000 



2000 



Bachiller et al. 



O' Linger et al. 



2003J). All are Class 



protostars except IRS 3A, which is a Class I protostar and bin a ry co mpanion to IRS 3B 
dCuriel et all Il99pl : iTerebev fc Padgettl Il997i lLoonev et all boool . l2003h . In the literature 
these sources have been identified with differing nomenclatures: L 1448- mm is also known as 
L1448C; L1448 IRS3 is also known as L1448N, and IRS 3A and 3B are also called N(A) and 
N(B), respectively. 

Outflows associated with star formation are a ubiquitous phenomenon. Observations 
and theory demonstr ate that the outflow originates within 100 AU of the central source 



[e.g. iShu et al.l Il995l ). With outflows originating deep within the circumstellar envelope , 
it is certain that the outflow will affect the evolution of the YSO (jArce fc Sargentl 120061 ). 
While the launching mechanism is still uncertain, the result of an outflow emanating from 
deep within the circumstellar envelope is the creation of a cavity in the envelope. As the 
protostar ages, the outflow cavity seems to widen. The widening is perhaps initially caused 
by precession of the outflow jet; howe ver the precession angles are found to be m uch narrower 
than the observed opening angle (e.g. Reipurth et al.|[2000 : Arce &: Sargent 2004 ) . The stellar 



wind must be the dominant mechanism for creating wide angle cavities (jArce &: Sargent 



20041 ) . As cavities are widened throughout the Class and Class I stages of star formation, 
the widening of the outflow cavities may act to dissipate t he envelope and halt mass infall 
(iShu et ali[l98~7l : IPadgett et al.lll999l : krce & Sargentl l2004h . 



We examine these outflow cavities because, in general, Class objects cannot be directly 
detected short ward of 10/im. The cavities allow light from the embedded central source to 
escape and scatter off dust in the cavity at near-infrared (NIR) wavelengths. NIR emission 
allows us to look into the cavity to the inner portions of the envelo pe. This bipolar cavity 



and outflow structure ha s been studied for more than 25 years (e.g lSnell et al.lll980l : iKaifu 
19871 : iMundt et al.l 119871 ): however, recent advances in instrumentation have enabl ed us to 
explore this structure in detail for Class sources (e.g. iNoriega-Crespo et al.l 120041 ) . In the 
past, only Class I sources could be studied in detail; high extinction in the circumstellar 
envelope greatly attenuates light at J, H, and Ks-bands for Class sources. The Infrared 
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Array Camera (IRAC) (IFazio et al.ll2004l ) on the Spitzer Space Telescope enables us to see 
through the dust enshrouding a Class source. For the first time, we are able to observe 
scattered light from Class sources with detail in multiple bandpasses. IRAC has channels 
numbered 1 through 4, corresponding to central wavelengths of 3.6/im, 4.5/xm, 5.8/xm, and 
8.0/xm wavelengths respectively. 



Radiative transfer codes (e.g. Whitney fc Hartmannll992 ; Lucas fc Roche 1997 ; Whitney et al 



200 3bf) can be used to model scattered light in the cayity from the embedded protostar (e.g. 



Lucas & Rochdll997l : IPadeett et al.lll999l : Eisner et alJl2005t iTerebev et al.ll2006l : IStark et al. 
20061 ) . at fluxes within the sensitivity of IRAC. Ground-based telescopes can detect scattered 



light emission in the J, H, and Ks bands with sufficient depth of field, requiring very long in- 



from more evolvec 


, Class I, sources (e.g. 


Whitnev & Hartmann 


1992. 


1993; 


Lucas & Roche 


1997; 


Eisner et al. 


2005; 


Stark et al. 


2006h. Class sources have not been studied in detail. 



The cavities as observed by IRAC will detect scattered light as well as thermal and molecular 
emission (if present), e.g. polycyclic aromatic hydrocarbons (PAHs), CO, and H2. Scattered 
light from the Class sources in L1448 is primarily only present in Ks and IRAC bands; J 
and H bands are very marginally detected. 

In this paper, we demonstrate that the idea of scattered light cavities is valid for Class 
sources in addition to Class I sources. Also, for the first time modeling of Class sources is 
performed primarily using the scattered light SEP an d images. Modeling these sources with 
the radiative transfer code by IWhitney et al.l (l2003bl ) enables us to determine geometric pa- 
rameters (inclination and opening angle) of the outflow cavity and loosely constrain physical 
parameters. To determine the true opening angle, we must simultaneously fit inclination. 
The parameters as constrained by the model are then compared to outflow observations 
at millimeter wavelengths. Finally, we discuss the idea of opening angle being a possible 
indicator of evolutionary state. 



2. Observations and data reduction 

2.1. IRAC and near-IR observations 

We observed the L1448 cloud with two separate pointings of IRAC on 2005 February 
25, one to cover L1448 IRS 3 and L1448-mm, and another to cover L1448 IRS 2. IRS 3 was 
observed in the High Dynamic Range mode with frame times of 30 seconds in 1 x 2 mapping 
mode and a cycled dither pattern of 31 positions and small scale factor achieving a total 
integration time of 930 seconds. The IRS 2 field was observed with the same parameters 
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but with 30 dither pattern cycle positions for a total integration time of 900 seconds. The 
difference in total integration times between the two fields is not significant. 

The L1448 fields had already been covered by the cores2disks (c2d) Cycle 1 Legacy 
program (jEvans et al.ll2003l ). but our science goals required much greater depth of imaging 
to be sensitive to scattered light in the outflow cavities. The additional depth has the added 
effect of a mitigating cosmic ray hits, producing a very clean final image with high signal to 
noise (S/N) ratio. This is necessary for the study of the outflow cavities, whose extended 
emission in these fields can be confused by the presence of background emission from the 
molecular cloud itself, as well as the HH objects associated with multiple YSOs in the region. 

Post-BCD pipeline products were solely used in this paper. Initially it was assumed that 
further reduction would be required, but after receiving the data, it was decided that the 
quality of IRAC data from the post-BCD pipeline was very high and sufficient for our science 
goals. Our data were processed by pipeline version S13.2.0. A color-composite image of the 
combined L1448 IRS 3 and IRS 2 fields is shown in Figure 1. Figure 2 shows a grayscale 
image of IRAC channel 2 with the positions of the major sources plotted. Close up views of 
the individual sources plotted with the PAs of CO outflows and source positions are shown 
in Figures 3, 4, and 5. 

Additionally, the paper makes use of data from t he Coordinated M olecular Probe Line 
Extinction Thermal Emission (COMPLETE) survey jRidge et alJboOfih . The NIR fields of 
L1448 were taken in J, H, and Ks bands, from the Calar Alto Observatory 3.5m telescope 
with the O mega2000 NIR prim e -focus imager. The observations and data reduction are 
desc ribed in Foster fc Goodman! (120061 ). a color-composite of the data are shown in Figure 
1 of iFoster fc Goodman! tood\ We calibrated the J, H, and Ks maps of L1448 to 2MASS 
standards (ISkrutskie et al.l 120061 ) . The error in flux calibration for J was 15%, H 11%, and 
12% for Ks. The inaccuracy of the pho tometric calibration may partly be due to the effect 
of cloudshine ([Foster fc Goodman! 120061 ). 



Photometric measurements of the fo ur sources were corrected for extinction using th e 
reddening law for infrared wavelengths (IRieke &: Lebofskyl Il985l ; llndebetouw et al.l 120051 ) . 
The extinction toward the sources comes from extinction maps of Perseus produced by the 
COMPLETE survey. The maps were made using 2MASS data and have 5' resolution. The 
extinction toward each source and their positions are listed in Table 1. 
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2.2. IRAC and near-IR photometry 

We performed aperture photometry on each source listed in Table 1 using the Im- 
age Reduction and Analysis Facility (IRAF). Photometry of the sources was not entirely 
straightforward since a background annulus could not be used in this situation due to ex- 
tended emission surrounding each source. Instead, using the imstat procedure in IRAF, we 
measured a large area of sky adjacent the source that was devoid of stellar emission due 
to high extinction to obtain an average background value per pixel. This value was then 
multiplied by the aperture area measured by IRAF and subtracted from the photometric 
value of our source. 

For each source, we measured a 1000 AU aperture radius centered on the positions listed 
in Table 1. Perseus is at a distance of roughly 250 pc, and the IRAC pixel scale is l'.'2/pixel, 
so we measured an aperture radius of 3.33 pixels. The pixel scale of the NIR images is 
O'.'45/pixel, making our aperture 8.89 pixels. We limited photometry to 1000 AU apertures 
to minimize contamination of the SED from extended H 2 emission, HH objects, and binary 
companions. Photometric values are listed in Table 2. 



3. Modeling 



To interpret our data, we used the radiative transfer code of IWhitney et al.l (l2003bl ). 
This is a multiple scatter and radiative equilibrium c ode using the Mont e Carlo method. 
The model assu mes a rotationally-fl attened envelope (ITerebey et al.lll984j ) with power law 
density r _3//2 (cf. lLooney et al.ll2003l ). We set a small centrifugal radius of 3 AU and an inner 
envelope radius of 3 times the dust destruction radius. A small centrifugal radius assumes 
little rotation and near spherical symmetry, and the dust destruction radius varies depending 
on particular central source parameters. The envelope also includes bipolar cavities; the 
cavity can be conical (streamline cavities) or they can be curved. We model curved cavities 
as this shape best fit our observations. Also, the model assumes a flared circumstellar disk, 
we use a disk with a radius of 50 AU; the inner radius is equal to the dust destruction 
radius. Embedded within the disk and envelope is the central source, defined by a stellar 
atmosphere and a given mass, radius, and temperature. The central source parameters may 
be somewhat unphysical for Class sources as the central source is still a protostar and not 
a pre-main sequence star. Observational constraints on the central source properties are not 
available for Class sources; therefore, we use the constraints that are available for Class 



I sources. iNisini et al.l (120051 ) measures effective temperatures of a few Class I sources to 
range from 3580 to 4900K and stellar radii of 1.2 to 2.9 R Q . These observed Class I central 
sources are similar to the central sources modeled in our grid. We assume a constant effective 
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temperature of 4000K and vary the stellar radius. 



With each run, the standard mode of the program outputs SEDs at 10 inclinations with 
bins centered at cos i = 0.05, 0.15, 0.95, with bin widths of 0.10. The code also has 



a "peeling-off" option (lYusef-Zadeh et al.lll984j ; IWood fc Reynolds! 1 19991 1 . which outputs a 
SED at only the specified inclination, without averaging over bins, and a high S/N image. 
The "peeling-off" option is much more computationally intensive than the standard mode, 
which is why we do not simply run the entire initial grid with this option. 

In modeling our data, we first compared our photometric measurements to SEDs from 
a grid of models run in the standard mode. Because of the large number of free variables, 
we only varied: luminosity, envelope mass, envelope radius, inclination angle, and opening 
angle, listed in Table 3. Opening angle refers to the half opening angle of the source; the 
angle as measured from the central axis to the cavity edge. To scale the stellar luminosity 
of the sources, we ran models with eight stellar radii (R*) at constant stellar temperature. 
Simply scaling a lower stellar luminosity model to higher luminosity does not give viable 
results because the inner disk and envelope radii depend on R*; scaling the luminosity 
would imply a higher temperature of the central source. However, accretion luminosity is 
inversely proportional to R*; we keep a constant accretion rate in our model grid rather than 
maintaining a constant accretion luminosity. 

This first grid was run in standard mode with only 1 million photons per run. This 
reduces computation time as there are 120 independent models per luminosity. Each model 
takes approximately 10 minutes to run on a 2.4 GHz Opteron processor with g77 compiled 
code. Once finished, we convolved the SED with filter functions for common bandpasses. We 
used a x 2 routine to obtain fits of the model SED to the observed SED at the IRAC and Ks 
bands taken with 1000 AU apertures; the sources are too weak at J and H bands. We do not 
include longer wavelengths because published photometric measurements are almost always 
upper limits due to their apertures encompassing multiple sources (our sources are all close 
binaries). Also, we are not attempting to fit disk properties with our grid as the envelope 
is most important in regards to scattered light. The disk properties are more important at 
l onger wavelengths e.g. submillimeter and mil l imeter, due to dust gr ain growth in the disk 
rtWolf-Chase et alil200ol : b'Alessio et~aD l200ll : Iwhitnev et alll2003bh . In Figures 6, 7, and 
8, we show the SED of the overall best fit model with our photometric measurements, we 
include photometric values from the literature in the SED for completeness. 

We define a good fit as being within the 90% confidence level. Our maximum allowed 
value is x 2 = 2, if we only fit IRAC bands, and 1.85 if Ks-band is included. We fit IRAC 
and Ks bands for L1448-mm and IRS 3B. IRS 2 was fit using only IRAC bands because it 
could not be fit at the 90% confidence level when Ks-band was included. Because IRS 2 is 
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the weakest source at Ks-band, it is possible that there is higher extinction towards IRS 2; 
this would affect Ks-band photometry more than IRAC photometry. IRS 3A was unable to 
be fit at the 90% confidence level with IRAC or Ks bands. These fits at the 90% confidence 
level are then taken to be inputs for detailed modeling. At the 90% confidence level, this first 
phase of modeling reduced our 9600 possible models to 36 models (see Table 4) for detailed 
modeling. The SED fitting does not tightly constrain physical parameters, however, we do 
observe certain trends from SED fitting. 

One trend that is observed from the SED fits at Ks and the IRAC bands is that lumi- 
nosity is not well constrained by our observed NIR SEDs. A model with the same physical 
parameters can appear with 3 or more luminosities. This occurs because Class YSOs 
emit the most flux in mid-infrared to submillimeter wavelengths. The luminosity is loosely 
constrained because fitting at only the shorter wavelengths is analogous to fitting a dog by 
its tail: you only roughly know if the dog is a great dane or chihuahua. Also, variations of 
opening and inclination angle can allow the SED to fit multiple luminosities. For example, 
a model that has a lower luminosity than the observed source can be a possible fit if the 
model has a wider opening angle and/or a smaller inclination angle. 

Another trend is the recurrence of inclination angle in multiple fits. For some sources, 
the same inclination angle is fit for several different luminosities, envelope radii, and masses. 
This reaffirms that the luminosity is not well constrained by our SEDs but also that the 
envelope radii and masses are not well constrained either. These trends are the reasons that 
this first phase of modeling only gives us a baseline for more detailed modeling. In order to 
confirm that a model fits the data, we must also compare a model image to the observed 
source. When the images are compared, many of the SEDs fit at the 90% confidence level 
are eliminated because images produced with the same parameters do not match the data. 

Our second phase of modeling used the "peeling-off" option of the Whitney code run 
with the parameter range of the SEDs fit at the 90% confidence level, producing a high S/N 
ratio image and an exact SED for a given inclination. The "peeling-off" option requires many 
more photons to produce a high S/N image. We ran our models with 10 million photons. 
Each model took approximately 6-8 hours on a 2.4 GHz Opteron processor. The images help 
us constrain the inclination angle and opening angle of the source. Different inclinations 
change the amount of flux measured from the source; flux increases from edg e-on to pole-on 



viewin g, and the appearance of the scattered light in the cavity changes (jWhitney et al. 



2003af ). 
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4. Discussion 

In modeling our sources, we had good success fitting SEDs and then refining these 
parameters with image comparison. However when we attempted to fit these images quanti- 
tatively to the data, we had mixed results. In general, the sources have complex morphologies 
when compared to our ideal models. Figure 9 displays images of the models which best match 
our data. The general trend is that the models overestimate flux in J, H, and Ks-bands, 
roughly match the flux in IRAC bands 1 and 2, and underestimate the flux in IRAC bands 
3 and 4. There are also morphological details of the data that cannot be fit by our models. 

Our models represent the most ideal scenarios of star formation, e.g. isolation, spherical 
symmetry, no external influences, etc. The data tell a different story. The sources we model 
have many complex morphological details not accounted for in an ideal scenario. All the 
sources are binary or multiple systems and the circumst ellar envelopes of th e sources are 



not necessarily spherical, the may be prolate or oblate (e.g lBonnell et al.lll996l ). In addition, 
molecular emission at NIR wavelengths is nearly always associated with YSOs in the form 
of HH objects, and emission from the disk or envelope around the protostar itself. Our 
models do not include flux from molecular lines present in NIR bands, PAHs, CO, and H2 in 
particular. H2 lines are present in all IRAC bands decreasing in intensity going toward longer 
wavelengths, the CO line is present in IRAC channel 2. On the other hand, PAH emission is 
not likely to be strong because the lo w mass YSOs in L1448 do not emit the far- u ltraviolet 



flux necessary to excite PAH emission (lAUamandola et al.lll989t iPeeters et al.ll2004l). Despite 



the p ossibility of line emission, spectra from the Infrared Spectrometer (IRS) (IHouck et al. 



20041 ). onboard the Spitzer Space Telescope, do not show evidence of line emission from any 
of the sources (IRS Disk Team, private communication). Unfortunately, IRS does not cover 
IRAC channel 2, thus we do not know if a CO line is present in any of these sources. 

Figures 1 to 5 show a large number of HH objects associated with outflows in close 
proximity to the sources. The presence of outflows and HH objects from several sources in a 
small area leads to the possibility of outflows influencing the circumstellar environment of a 
source. Because our sources do not exist in the ideal environment of the models, we are not 
able to quantitatively compare their images. Nevertheless, despite the many ways in which 
the sources do not fit an ideal scenario, the models do generally describe the circumstellar 
structure of the protostars in L1448, and a qualitative image analysis can be performed. 

Our primary method of matching model images to the data is the comparison of flux 
versus position angle at a given distance from the central source. We use an annulus of 
radius 3000 AU (10 pixels) from the source and a thickness of 600 AU (2 pixels). The flux 
density is measured at 64 discrete locations, or segments, along the annulus by a custom 
FORTRAN program. Only full pixels within a segment are used to measure the average 
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flux density per segment, not partial pixels, so there is some flux lost in this measurement 
process. However, the loss is not crucial since the same loss is in both the models and the 
data, and the plots are only used for qualitative comparisons rather than numerical fitting. 

For comparison purposes, we rotate all of our cavity images such that their outflow 
position axis is vertical, 1.57 radians (90°). Starting from a position of zero radians at the 
horizontal, the top of the image is positive angle and the bottom is negative angle. The 
center of the cavities are located at approximately 1.57 radians for the blue-shifted side and 
-1.57 radians (-90°) for the red-shifted side. Depending on the presence of confusing sources, 
we may not have used the full angular range for fitting. From the plots shown in Figures 10 
to 12, it is clear that our model flux values do not entirely match the data flux values. 

The most important features we examined for model and data agreement were: how 
well the boundaries of the cavity match in the models and data and how well the amount 
of flux agreed in the IRAC bands for these plots. The boundaries of the cavity are where 
the flux starts to continually increase from the zero level up to the peak and then drop back 
to the zero level. A model with the correct combination of opening angle and inclination 
will match the boundaries in the data. A larger or smaller opening and/or inclination angle 
will not match the data and not be considered a fit. There is some degeneracy between the 
inclination and opening angles, however, we use the SED fitting to disentangle them. 

We then looked at how well the flux from the cavity at its peak levels matched the 
models. The general pattern of the flux in the models is to rise to a peak level, remain 
approximately at this level for the extent of the cavity, and then decrease back to zero. This 
pattern is observed with IRS 3B, but IRS 2 and L1448-mm rise to a peak in the center of the 
cavity, then drop back, covering the same angular extent as the model. The models do not 
simulate this feature. It may be evidence of more scattering at the center of the cavities of 
IRS 2 and L1448-mm than at the cavity-envelope interface, this is discussed in more detail 
below. For L1448-mm and IRS 2, we looked at how well the flux at the peak overlapped 
the maximum flux level of the model. This method does not quantitatively constrain the 
physical parameters, but the SED does quantitatively fit the physical parameters. Matching 
the images is simply a further constraint on the SED fit. It is fair to say that there could 
be more possible SED and image fits within the resolution of the model grid; we adopt the 
resolution of the model grid as our error bars for inclination and opening angle. Table 4 
describes the possible uncertainty of fit parameters in more detail. 

As mentioned previously and shown in Figures 10 to 12, there are inconsistencies be- 
tween the model and data that are dependent upon wavelength. A possible reason for these 
inconsistencies is that there may be more extinction toward the sources than derived from 
published extinction maps. We base this assertion on the fact the models over produce flux 
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at the J, H, and Ks-bands. We tried different variations of Ay attempting to match the 
J, H, and Ks-bands. The appearance of the models began to better agree with the data if 
we apply an Ay = 20. However as we better match the NIR bands, our models and data 
diverge in the IRAC bands. Also, at high extinctions, our observed SEDs do not fit any 
models. An alternative explanation for these inconsistencies could be incorrectly assumed 
dust properties. 

Incorrectly assum ed dust properties seem more likely than an extinction of Ay > 20. 



Whitney et al.l (j2003bl ) uses interstellar medium (ISM) dust models in the cavity and enve- 



lope. An ISM dust model may not be completely reali s tic for protostellar envelopes. The 



dust albedo of the ISM models used by IWhitney et al.l (l2003bl ) increases steadily at wave- 



lengths short ward of 10/im and peaks around 1/im. If this and the published values of 
Ay are correct, there should be much more scattered light than we observe in the J and H 
bands, and less scattered light in IRAC channels 3 and 4. We do not observe this trend. 
Also, we do not know the size of the dust grains in the observed cavities or the effects of dust 
grains varying in size. Larger dust grains could cause more scattering in IRAC channels 3 
and 4 and less scattering in the J and H bands. Additionally, cavity density is assumed as a 
constant in our models and not varied in the grid, however we did explore the effects of half 
an order of magnitude and an order of magnitude higher or lower cavity density than listed 
in Table 3. There is a very important caveat to modeling a constant cavity density, partly 
depending on the opening angle. A cavity density an order of magnitude greater than our 
assumed value begins to have a mass that is not negligible as compared to the mass of the 
envelope. 

A greater cavity density than our assumed value (see Table 3) decreases scattered light 
at short wavelengths (J, H, Ks, and IRAC channels 1 and 2). A lower cavity density than 
our assumed value does not increase or decrease scattered light at short wavelengths. IRAC 
channels 3 and 4 are not affected by higher or lower cavity density; the model predicts that 
nearly all emission in these two channels is thermal due to the ISM dust model used for the 
cavity and envelope. From examining the effects of varying cavity density, the models are 
predicting that most scattered flux should come from the envelope-cavity interface rather 
than dust grains in the cavity; the dust grains in the cavity seem to extinct flux more 
than they scatter it. As discussed previously, the model does not necessarily reflect the 
morphology of scattered light in our data. Our SED fitting at a 1000 AU aperture is not 
greatly affected by higher or lower cavity density. However, an SED measured with a larger 
aperture (e.g. 5000 AU) and a lower cavity density will measure more J, H, and Ks-band 
flux, conversely higher cavity density will emit less J, H, and Ks-band flux due to extinction. 
In addition to possible variations in cavity density, the cavity density is not likely constant 
as we assume. It could be decreasing as a function of distance from the central source, or 
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vary as function of an gle within the cavity to reflect high and low velocity outflow regions 
( IWhitney et al.l l2003bl ) . Limb brightened cavities are det ected in low velocity CO emission 
and the central jet is detected as high velocity emission (jBachiller et al.lll995l ). We do not 
explore these possible cavity density variations. 

Despite some of the obvious limitations of the images only qualitatively describing our 
observations, the physical parameters derived from modeling are realistic and conform to 
previ ous work (e.g lBachiller et aLlllQQSKlO'Linger et al.lll999l ; lLooney et al.ll2003l ; lKwon et al. 



20061 ). In areas where our results do not conform, there are reasons and new insights for 
the discrepancies highlighted by our data. A caveat that we must mention is that the best 
fit envelope radius we quote may not be a constraint in reality due to the SED only being 
taken at a 1000 AU aperture. We regard the distance the scattered light extends as a 
better measurement of envelope radius, though this value is probably a lower limit because 
cavity density may vary inversely with radius. The most realistic constraint of envelope 
radius would be a high resolution extinction map derived from IRAC data; however, such an 
analysis this is beyond the scope of this paper. For all the sources, except IRS3A, we were 
able to match a model to the data using our technique of SED fitting and image matching 
(Table 4). In the following subsections, we discuss each individually. 



4.1. L1448 IRS 2 



In many ways, L1448 IRS 2 is the most ideal source in the field. It has a curved outflow 
cavity in which both the red and blue-shifted sides are clearly detected (see Figures 3 and 
9a). Also, at first glance, it appears to have evolved in relative isolation compared to the 



inter acting outfl ows of IRS3A, I RS 3B, and L1448-mm (IKwon et al.l 120061 ; lO'Linger et al. 



20061 ). Although IVolgenau! (120041 ) detect IRS 2 as a compact binary source, our IRAC data 
and the MIPS data from the c2d project, do no t detect any obvious in dicator of its binary 
nature. The companion lies on the PA (~ 138°) (IWolf-Chase et al.ll2000l ) of the outflow from 
IRS 2 separated by ~2500 AU (~ 10"). We will refer to the primary source as IRS 2A and 
the companion as IRS 2B. We are still able to qualitatively fit an image to IRS 2A as IRS 
2B does not appear to affect the morphology of the IRS 2A in our data. 

As shown in Figures 9a and 10, the modeled image approximately defines the same 
cavity boundaries as our data but quantitative differences are evident. First, there is more 
scattered light emission in the model than the data at 3.6/zm and only one cavity boundary 
is well defined in the data. The other boundary is indeterminate at 3.6 /im. At wavelengths 
longer than 3.6 /zm, the flux in the middle of the cavity dominates the features. However, 
at 5.8 and 8.0/mi the red shifted side of the cavity shows limb brightening. Also, we observe 
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the cavity flux increasing with wavelength, while the models predict that the cavity flux 
will fall at longer wavelengths, due to properties of the dust grain model. As discussed 
previously, incorrectly assumed dust properties may account for the wavelength dependent 
discrepancies. Also, H 2 and PAH emission are not likely to be present in the cavity as IRS 
spectra do not show evidence of emission lines (IRS Disk Team, private communication). It 
is also important to note that the emission peaks in each band fall approximately on the PA 
of the outflow from IRS 2, and the annulus used to construct Figure 10 is coincident with 
the location of IRS 2B. 

Based on the compared images and SED fit for IRS 2A, the inclination and opening 
angles are 57°±g and 20° ±5 respectively. Our data show that scattered light extends ~ 7000 
AU on the blue-shifted side and ~ 5750 AU on the red-shifted side. The model image shows 
scattered light extending to 5000 AU on the red and blue-shifted sides of the cavity. It is 
unclear whether or not the scattered light will end at the edge of the envelope as it does in 
the models. 

As stated above, IRS 2 is a compact binary; our data only resolve the cavity of IRS 2A. 
Our IRAC observations could have detected another cavity associated with IRS 2B, which 
would show the direction of its outflow. The outflow direction of IRS 2B is impor t ant as 
uncertainty remains as to the source of the eastern outflow, HH195E. iBally et al.l (119971 ) 
infers that the source of HH195E is L1448 IRS 1, a Class I source, because o f an apparent 
east-west outflow centered at its location. However, I Wolf- Chase et al.l (120001 ) suggest that 
another undetected source might be responsible for this outflow. The lack of an observed 
cavity associated with IRS 2B, does not exclude an outflow source from being present. There 
is a Class source in th e region for which we do not clearly detect a c avity, L1448 NW also 
named L1448 IRS 3C JTerebev fc Padgettl fl~99~7l : lLoonev et all boooh . It seems probable 
that only high resolution mapping of molecular outflow tracers will fully disentangle these 
outflows. 



4.2. L1448 IRS 3B 



As Figures 1 and 4 show, IRS 3B is a confused source, lying in close proximity (~ 7") 
to IRS 3A, its binary companion. In additio n to being a confu sed source, the cavity is not 
symmetric about the outflow PA (~ 105°) (IKwon et al.ll2006l ). In an ideal situation, the 
cavity would be symmetric about the PA of molecular outflows. With this PA as the axis 
of symmetry, the southern angle of the cavity is greater than the northern angle, see Figure 
4. Only one side of the cavity is clearly visible in our data, making the asymmetry and 
PA difficult to confirm. However, our data support a ~ 105° PA because HH objects are 
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detected southeast and northwest of IRS 3B on this PA. 



Due to this asymmetry, no model fills the complete extent covered by the observed 
cavity; models match the northern edge of the cavity, but not the southern edge. One 
possible reason for this asymmetry is that IRS 3B is being affected by the outflow of L1448- 
mm. As seen in Figure 1, the outflow jet from L1448-mm is curving toward IRS 3B, possibly 
confusing the outflow or perhaps eroding that side of the envelope causing the cavity to widen. 
O'Linger et al.l (120061 ) posit that the L1448-mm outflow eroded the envelope of IRS 3 A. It 
may be possible that this is ha ppening to IRS 3B as well. Another possibility is that the 
outflow of IRS 3B is precessing (jGueth et al.lll996l ; IShepherd et al.l I200CH : lArce fc Goodman 
200ll ; iLebron et all l20od : lYbarra et aD l2006h . HH objects are observed on the PA of IRS 
3B and others are present off this PA. If these HH objects located off the PA are associated 
with IRS 3B, a precessing outflow could account for the HH objects on two PAs as well 
as an asymmetric cavity. A precessing outflow is not something our data can confirm; it 
merely suggests the possibility. We must note that L1448 NW (IRS 3C), a Class source, is 
present near IRS 3B. Despite remaining undetected in our data, it appears to be driving an 
outflow on a PA of ~ 128°. This PA is estimated from tracing HH objects which appear to 
be collimated toward L1448 NW; molecular line observations have not detected an outflow 
from L1448 NW. Incidentally, there are HH objects present in the c2d data roughly ~ 0.7 
pc distant from L1448 NW which appear to be roughly coincident with the PA of ~ 128°. 
However, the large distance makes a definite source association uncertain. Nevertheless, it 
is clear that possibility of multiple outflow sources confuses our analysis of the outflows in 
association with HH objects. 

As shown in Figures 9b and 11, the model predicts flux from the red-shifted side of the 
cavity at levels too low to distinguish. HH objects with strong emission are present in the 
vicinity of where the red-shifted side would appear. The blue-shifted side is relatively well 
matched by model predictions in IRAC channels 1, 2, and 3. Despite being asymmetric, 
the flux emitted by cavity of IRS 3B corresponds better to its best fitting model than the 
other sources. IRS 3B is also the only source that does not peak in flux at the center of the 
cavity in any band. Figure 11 shows some evidence of limb brightening in all bands. Even 
though IRAC matches very well, discrepancies at H and Ks bands are still present. These 
discrepancies may be explained by incorrect dust properties and cavity density as described 
in §4. H2 and PAH emission are not likely to be present in the cavity as IRS spectra do not 
show evidence of emission lines (IRS Disk Team, private communication). 

Even though IRS 3B is not an ideal source to model, as we do not observe both sides 
of the cavity; the parameters of the model which best describes the data are consistent with 
previous studies. The inclination and opening angles, 63°±g and 20° ± 5 respectively, agree 
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with parameters derived in lKwon et al.l (120061 ). In the IRAC data, scat tered light extends to 
~ 7500 AU, possibly implying an envelope radius of the same extent. lLooney et al.l (120031 ) 
derived an envelope radius of 8000 AU, however, the adopted distance was 320 pc. Correcting 
this value for our distance (250 pc), the observed envelope radius is 6250 AU. It is important 
to note that this envelope is shared with IRS 3A, however, each source appears discretely in 
our data and can be analyzed independently. 



4.3. L1448-mm 

L1448-mm also has a complex morphology that makes modeling somewhat more diffi- 
cult. In Figure 5, you see that the source has two bright spots of emission. One is directly 
north, associated with its cavity and the other is directly south. The central source is located 
adjacent to the northern brightness peak. The separation between the northern and southern 
peak is ~ 7" (~ 1800 AU) in IRAC and MIPS channel 1, see Figure 13. Only the northern 
peak is observed in K s-band. These data suggest that L1448-mm may be binary in nature. 



Bachiller et al.l (119951 ) mapped L1448-mm at A = 2.7mm continuum with a resolution of 3 



the data s how t hat L1448-mm is extended along the molecular outflow PA (~ 1 59°), and 



Volgenaul (120041 ) detect three distinct sources at A = 2.7mm. The conclusion of I Volgenau 
( 120041 ) was that L1448-mm is a multiple system. Our data support this conclusion, detecting 
one companion source in our IRAC and MIPS data. We will refer to the northern source, 
as L1448-mm A and the southern source as L1448-mm B. Only L1448-mm A is associated 
with the scattered light cavity. 

Despite this complex morphology, the scattered light image from L1448-mm A is ap- 
proximately matched by the model, see Figures 9c and 12. The model matches the data 
best in IRAC channel 1. In IRAC channels 2, 3, and 4, the flux of the model is too low as 
compared to the data. However, aside from the indentation on the west side of L1448-mm 
A, the models approximately match the cavity boundaries of the data in each channel. The 
data in Figure 12 may show evidence of more scattering at the cavity center, as the flux 
approximately matches the boundaries predicted by the models but peaks in the middle of 
the cavity. Notice that the peak flux also lies on the PA of the outflow. Also, some of the 
wavelength dependent inconsistencies could be the result of improperly assumed dust prop- 
erties and/or cavity density as discussed previously. No obvious emission lines are present 
in this source in the wavelengths covered by IRS (IRS Disk Team, private communication). 

Although the models do describe the general morphology, the models were not able 
to completely match the details of L1448-mm A. The northern side of the cavity is ap- 
proximately matched by the model; however, there is no southern cavity in our data, only 
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L1448-mm B. The brightness of source B could be masking the southern side of the cavity. 
Also, the scattered light in the model does not extend as far as in the data. The data appear 
to show scattered light out to ~ 9500 to 8500 AU from source A. There is a jet present in this 
region and the emission extending out to 9500 AU could be light scattering on the ambient 
mediu m. Also, the limb brightened cavity walls observed in CO J = 1 — > by lBachiller et al. 
( 119951 ) extend ~ 8800 AU on the northern side and ~ 10000 AU on the southern side. 



The inclination and opening angles of the best matched model, 49° ± 8 and 15° ± 5 
respectively, are different, but not extraordinarily so, from the inclinat ion angle of ~ 70° an d 
opening angle o f 22.5° derived from molecular outflow observations (jBachiller et al.lll995l ). 
Girart fc Acordl (120011 ) derives nearly the same r esult with a dif f erent m ethod. A cavity 
i mage using the opening angle and inclinations from lBachiller et al.l (119951 ) or lGirart fc Acord 
(120011 ) could approximately match the data; however, a model SED with these parameters 
was not fit at the 90% confidence level. 

The fact th at L1448-mm is a multiple system composed of at least 3 sources, possibly 4 
(jVolgenaull2004l ) . complicates the derivation of physical parameters. Because the companion 
source is so close in proximity to L1448-mm, our SED measurements are slightly contami- 
nated with flux from L1448-mm B, adding uncertainty to our results. In §3 it is mentioned 
that flux increases with as the inclination gets smaller. The additio nal flux from source B 
sour ce is likely the reason for fitting smaller inclination angles than iBachiller et al.l (119951 ) 
and iGirart fc Acordl (bpOlh . 



4.4. L1448 IRS 3A 



Despite our success with the other sources, L1448 IRS 3A was not constrained at the 
90% confidence level by our SED fitting. This however, is not unexpected because IRS 3A is a 
Class I source with an envelope mass of ~ 0.09 — 0.29M ? (lLooney et al.ll2000i ; lO'Linger et al. 



20061 ). Our model grid does not attempt to model the Class I stage envelope masses or 
opening angles. Additionally, the disk properties of a Class I source should have different 
model input parameters than a Class source. A C lass I source is dominated by the disk 
more than the remnant envelope (ILooney et al.ll2000i ). 



Despite not being constrained by modeling, some general statements can be made from 
the data. The peak emission in IRAC channel 1 is very near the centr al source location . 
Because the source is thought to lie at an inclination of approximately 90° (IKwon et al.ll2006l ). 
an emission peak within a few hundred AU from the central so urce could be evidence of an 
envelope that is nearly dissipated and a wide opening angle. lO 'Linger et al.l (120061 ) posit 
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that the envelope may have been eroded by the outflow from L1448-mm. This hypothesis 
is supported by the apparent lack of scattered light on the southern side of the source, see 
Figure 4. If the envelope has been eroded away on the side facing L1448-mm, there would 
be no scattered light on t hat side of the source because the cavity in the envelope is the 
primary scattering surface (jWhitney et al.ll2003bl ). 

Additionally, we can look for HH objects on the PA of the molecular outflows. In the 
case of IRS3A, we find emission from HH objects to the north and south of the source that 
lie on the PA (~ 155°) (IKwon et al.ll2006l ). see Figure 1. The HH objects to the south appear 
to be collimated toward IRS 3A with conical symmetry. This is probably not scattered light 
emission because the appearance of scattered light in our data is smooth and continuous 
while these objects appear "lumpy." The HH objects that appear associated with IRS 3A 
are all within 50". 



4.5. Opening angle and evolutionary state 



We stated in §1 t hat previous work posit that the cavi t y opening angle will increas e 



Class I protos t ars m odeled by Eisner et al 



as the protostar ages flShu et all Il987l: IPadgett et all 1199a lArce fc Sargentl 12004 120061 ). 



(bo05h have opening angles of 25 - 28° and 



Terebey et all (120061 ) constrains the opening angle of TMC-1, a Class I protostar, to be 
40°±5. Additionally, a small Hub ble Space Tele s cope (HST) survey of near edge-on Class 
I objects in Taurus performed by IPadgett et al.l (119991 ) measures apparent opening angles 
> 30° for all but one source. However, in measuring these opening angles the authors do 
not account for inclination effects which could account for the outlier. In their sampl e , they 
find that opening angle is inversely proportional to circumstellar mass. IStark et al.l (120061 ) 
model these sources in Taurus and obtain fits for opening angles of 30° (DG Tau B), 25° 
(IRAS 4248+2612, 04301+2247), and 20° (IRAS 04016+2610, CoKu Tau/1). 

Our models of the Class sources in L1448 conform to opening angles of 20° and 15°. If 
the opening angle does indeed widen with age, there may be a small range of opening angles 
for which the transition from Class to Class I occurs. With the data to date, we suggest 
that this transition range is probably between 20° to 30°. This result, inferred from infrared 
ima ging of cavities, is con si stent with studies of o utflow cavities at millimeter wavelengths 
(e.g lArce fc Sargentl 120061 ) . lArce fc Sargentl (120061 ) observes that opening angles of Class 



source s to be less than 27?5, and Class I opening angles are wider. Also, lArce &: Sargent 



(120061 ) propose an empirical model for the outflow-envelope interaction which causes the 
opening angle to widen with age. Though our dataset small, it sets the ground work for a 
more comprehensive study of opening angles observed from cavities in IRAC wavelengths. 
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5. Conclusions 

We have demonstrated that the general rule of thumb "Class sources are invisible 
at A < 10/im" is incorrect. New advances in instrumentation allow us to probe into the 
circumstellar envelope by imaging the scattered light from the central protostar that is 
escaping out the cavity, carved by the powerful molecular outflows. In fact, the scattered 
light emission is brightest in the IRAC bands due to the lower extinction. These images 
impose, for the first time, constraints on the central protostar and environment in a new and 
exciting way. In particular, we have shown that the general morphologies of Class objects 
(i.e. circumstellar envelope with a cavity, circumstellar disk, and embedded protostar) are 
consistent with our scattered light images. The ability to which the models match the data 
depends largely on the circumstellar environment of the protostar. If the protostar has a 
very ideal morphology (i.e. spherical and axial symmetry), it is possible to well model the 
source. Complex morphologies e.g. IRS 3B and L1448-mm are more problematic. Even so, 
an ideal model that approximates the structure of an non-ideal source provides important 
constraints on its circumstellar structure, mainly opening angle and inclination angle, as well 
as providing loose constraints on the source parameters. 

We are able to fit an opening angle and inclination angle to our sources using both 
images and SED fits. For L1448 IRS 2, the source that best exhibits a bipolar cometary 
nebula, our best fit inclination angle is 57°±g and the opening angle is 20°±5. L1448 IRS 3B, 
a source which may have an asymmetric cavity, has a best fit inclination angle of 63°±g and 
the opening angle is 20° ± 5. L1448-mm, the source that drives one of the most prominent 
examples of molecular outflows, has a best fit inclination angle of 49° ± 8 and the opening 
angle is 15° ± 5. 

Deriving the outflow opening and inclination angles are crucial to understand the 
sources. There is a probable link between the age of the source and opening angle; ob- 
servations and models suggest that the opening angle widens with age. The inclination 
angle can completely change the appearance of the cavity, interacting with the apparent 
opening angle. In other words, one must derive the true opening angle of an embedded 
source with the inclination to decouple the interdependency of the geometry and extinction. 
We assert that the opening angle of a protostar can be indicative of its evolutionary state 
based on our modeling of Class sources in L1448 and previous studies of outflow cavities. 
Also, comparing the inclination angle derived from modeling to an inclination angle derived 
from molecular outflow observations is an important baseline. A possible constraint on en- 
velope radius could be the distance the scattered light extends from the source; however, 
this distance is likely a lower limit. 

Spitzer gives us a very unique and detailed glimpse of circumstellar structure around 
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protostars. At the same time, however, it reveals that the circumstellar structure and for- 
mation environments are generally far from being idealized cases. Analysis of models versus 
data does provide important constraints of the circumstellar structure of protostars, but we 
may still be missing important details. In the near future, we may be able to more quan- 
titatively c onstrain the circumstellar structure. New t hree-dimensional radiative transfer 
codes (e.g. IWhitney et al.ll2005l ; llndebetouw et al.ll2006l ) may enable us to model even the 
most unique circumstellar environments. However, with new radiative transfer codes the 
multitude of parameters increases as does the time required to model a comprehensive grid. 
Observational programs must be able to constrain input parameters to reduce complexity. 
Future instruments, e.g. James Webb Space Telescope, Atacama Large Millimeter Array, will 
be able to provide astronomers with highly sensitive and detailed observations of protostars, 
building upon the foundation that Spitzer has laid. 
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Table 1. Source List and Extinction 



Source Coordinates (J2000.0) 
Right Ascension 



Declination 



A a 



L1448 IRS 3 A 
L1448 IRS 3B 
L1448-mm 
L1448 IRS 2 



03h 25m 36.532s 
03h 25m 36.339s 
03h 25m 38.8s 
03h 25m 22.4s 



+30°45m 21.35s 4.57 

+30°45m 14.94s 4.57 

+30°44m 05s 5.4 

+30°45m 12s 4.37 



a Derived from lRidge et al.l (120061 ) . 



Table 2. Photometry 



Source 


Aperture 


Fj 


Ftf 


F Ks 


F3.6 


F4.5 


F5.8 


Fs.o 




(AU) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


L1448 IRS 3 A 


1000 


0.172 


0.293 


0.750 


2.54 


16.6 


51.6 


96.6 


L1448 IRS 3B 


1000 






0.380 


2.44 


6.85 


9.80 


11.7 


L1448-mm A 


1000 




0.144 


1.07 


4.80 


15.5 


29.7 


65.1 


L1448 IRS 2A 


1000 






0.117 


1.15 


5.38 


10.2 


15.8 



Note. - Photometric measurements have been corrected for extinction. J, H, and Ks, 
flux values have uncertainties of 15%, 11%, and 12% respectively. IRAC flux values all have 
uncertainties of 15%. 
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Table 3. Model grid parameters 



Parameter 


Description 


Value 






9 HQ ^ ^ R9 ZL 1 8 A P.7 ^19 ^ ^ ^ Q1 
z.uy, o, o.oz, i.io 'i.u / , o.iz, o.oo, o.yi 






zinnn 


T (J \ 


aystem luminosity 


9 11 9 8/1 fi^ /I W K « AO 7 AK O AO 
Z.LL, Z.041, O.OO, 41.01, 0.01, O.^to, (.410, O.^tZ 


M (M } 


OLeildl lildSS 


u.o 


A/r /A/r \ 


Disk mass 


U.Uo 


a 


Disk radial density exponent 


2.25 




Disk scale height exponent 


1.25 


M disfe (M yr- 1 ) 


Disk accretion rate 


2 xl(T 7 


Retesfc,mm(Rs«f)) 


Disk inner radius 


1 


R-disfc,maa; ( ATJ) 


Disk outer radius 


50 


R C (AU) 


Centrifugal radius 


3 


Rem;, ram (Rsub) 


Envelope inner radius 


3 


Renu , mai ( AU ) 


Envelope outer radius 


5000, 6500, 8000 


M c 


Envelope mass 


0.25, 0.5, 0.75, 1.0, 1.25, 1.5, 1.75, 2.0 


@open 


Opening angle 


10°, 15°, 20°, 25°, 30° 


@inc 


Inclination angle 


cos i = 0.05, 0.15, 0.95 


p c (g cm' 3 ) 


Cavity density 


3.34 x 10~ 19 



a The system luminosity is not an explicitly varied parameter, but it is greatly dependent on 
R*. 

h R su b is the dust destruction (sublimation) radius. This is calculated internally by the model 
and will vary with the stellar radius/luminosity. 

c The envelope mass is set by the mass infall rate. The mass infall rate for a given envelope 
mass changes with different envelope radii. For brevity, we only list the approximate envelope 
masses. 
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Table 4. Best Fit Parameters 



Source 



Luminosity 



(L ) 



Stellar b 
Radius 

(Ro) 



Envelope c Envelope d Opening 
Radius Mass Angle 

(AU) (M ) n 



Inclination 
Angle 

(°) 



SED 

Xr 



L1448 IRS2 4.6 
L1448 IRS 3B 4.6 
L1448-mm 7.5 



4.2 
4.2 
5.5 



5000±1500 
5000±1500 
6500±1500 



1.50±0.5 
1.00±0.5 
1.75±0.5 



20±5 
20±5 
15±5 



57±| 
63± 7 6 
49±8 



0.374 
0.603 
1.62 



a The value found to best fit the 1.66 - 8.0/xm data; it is likely a lower limit on the bolometric 
luminosity. 

b This value depends on assumptions in our model grid and best-fit luminosity. 

c Envelope radius is weakly constrained by the fits. The error bars represent fit ranges within the 
model grid which likely underrepresent the total uncertainty. 

d The uncertainty in envelope mass is tied to the u ncertainty in envelope radius. Also if the density 
power law is different from r~ 3 / 2 ( Loonev et al. 2003 ) the envelope mass is more uncertain. A fit within 
the model grid simply fits the optical depth necessary to produce a similar SED for a given luminosity. 
As seen in Table 5, a multitude of different envelope masses and radii produce an SED fit. 
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Table 5. Results from SED fitting 



Source Luminosity Stellar Envelope Envelope Opening Inclination SED 

Radius Radius Mass Angle Angle Xr 

(L ) (R ) (AU) (M ) (°) (°) 



L1448 IRS 3B 



L1448-mm 



L1448 IRS2 



2.11 


2.09 


8000 


1.25 


15 


63 


1.27 


2.84 


3 


8000 


1.75 


20 


63 


1.42 


4.57 


4.18 


5000 


.50 


10 


81 


1.55 


4.57 


4.18 


5000 


1 


20 


63 


0.603 


5.51 


4.67 


8000 


1.25 


10 


81 


1.59 


6.48 


5.12 


5000 


1.75 


20 


57 


1.19 


6.48 


5.12 


6500 


1 


10 


76 


0.555 


6.48 


5.12 


6500 


1.5 


10 


49 


1.29 


6.48 


5.12 


6500 


1.75 


20 


70 


1.84 


6.48 


5.12 


8000 


1.5 


10 


63 


1.41 


6.48 


5.12 


8000 


1.75 


15 


81 


1.95 


6.48 


5.12 


8000 


2 


15 


70 


1.86 


7.45 


5.53 


5000 


1 


15 


81 


1.53 


7.45 


5.53 


6500 


2 


15 


57 


0.909 


8.42 


5.91 


8000 


2 


15 


81 


1.77 


2.11 


2.09 


5000 


.75 


10 


32 


1.70 


2.11 


2.09 


6500 


.75 


15 


57 


1.20 


2.11 


2.09 


6500 


1.25 


30 


63 


1.79 


2.11 


2.09 


8000 


1.75 


30 


63 


0.822 


2.84 


3 


5000 


2 


15 


32 


1.85 


2.84 


3 


8000 


1.25 


10 


41 


1.97 


3.65 


3.62 


6500 


1 


10 


41 


1.91 


3.65 


3.62 


6500 


1 


15 


57 


1.70 


4.57 


4.18 


6500 


1.5 


20 


57 


1.63 


4.57 


4.18 


8000 


1.25 


10 


49 


1.64 


5.51 


4.67 


5000 


.75 


10 


41 


1.72 


5.51 


4.67 


5000 


1 


15 


49 


1.82 


5.51 


4.67 


8000 


1.75 


10 


41 


1.95 


7.45 


5.53 


6500 


1.75 


15 


49 


1.62 


8.42 


5.91 


5000 


1.75 


25 


57 


1.21 


2.11 


2.09 


5000 


1.75 


20 


49 


1.97 


2.11 


2.09 


6500 


2 


30 


63 


1.93 


2.11 


2.09 


8000 


1.75 


15 


57 


1.03 
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Table 5 — Continued 



Source Luminosity Stellar Envelope Envelope Opening Inclination SED 
Radius Radius Mass Angle Angle Xr 

(L ) (R ) (AU) (M ) (") (°) 



3.65 


3.62 


5000 


1.25 


15 


49 


0.752 


4.57 


4.18 


5000 


1.5 


20 


57 


0.373 


6.48 


5.12 


6500 


2 


10 


41 


1.79 
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Fig. 1. — IRAC image of L1448, color composite of Channels 1 (blue), 2 (green) and 4 (red). 
This image is combined with the c2d data for full wavelength coverage across the field of 
view. 
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Fig. 2. — IRAC channel 2 image of L1448. Major sources are marked: IRS 2 is marked with 
an open box, IRS 3 is marked with a X, and L1448-mm is marked with an open diamond. 
The companion sources are identified in Figures 3, 4 and 5. 
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Fig. 3. — Zoomed in images of cavities for IRS 2. Images are 2' (30,000 AU) on each side. 
The solid line is the PA of the CO outflow from IRS 2A from I Wolf- Chase et all feoooh . the 
cross is the location of the peak millimeter emission from IRS 2A, and the star is the location 
of IRS 2B. Note that there is no obvious detection of IRS 2B. 
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Fig. 4. — Zoomed in images of cavities for IRS 3B. Images a re 2' 30,000 AU) o n each side. 
The solid line is the PA of the CO outflow from IRS 3B from lKwon et al.l (120061 ). The cross 
is the location of the peak millimeter emission from IRS 3B, the asterisk is the location of 
IRS 3A, and the star is the location of IRS 3C (L1448 NW). Note that the cavity of IRS 
3B is not symmetric about the PA of the CO outflow. Also, there is a marginal detection of 
emission near IRS 3C in IRAC channels 3 and 4. 
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Fig. 5. — Zoomed in images of cavities for L1448-mm. Images are 2' (30,0 00 AU) on each 



side. The solid line is the PA of the CO outflow from L1448-mm A from Bachiller et al. 



( 119951 ) . The cross is the location of the peak millimeter emission from L1448-mm A, and the 
star is the location of the peak millimeter emission from L1448-mm B. Note that L1448-mm 
B is undetected in the Ks-band and becomes more prominent at longer wavelengths. 
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Log A(/xm) 



Fig. 6. — SED of of best fit model for IRS 2A at apertures of 1000 (dashed line) and 5000 
AU (solid line), plotted with photometry. Data points represented as boxes were included 
in fitting, triangles were not used in fitting but are included for completeness. Wavelengths 
greater than 8.0/im are upper limits, as the aperture sizes are greater than our 1000 AU 
aperture. IRS 2 is a multiple system and all data except the A = 2.7mm data point include 
the companion sourc es. A = 12, 25, 60, 100, 450, 850, and 1300/mi d a ta po ints are from 



O 'Linger et al.l fjl999f ) . and the A = 2.7mm data point is from IVolgenaul (12004 ). 
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Fig. 7. — SED of of best fit model for IRS 3B at apertures of 1000 (dashed line) and 5000 
AU (solid line), plotted with photometry. Data points represented as boxes were included 
in fitting, triangles were not used in fitting but are included for completeness. Wavelengths 
greater than 8.0/im are upper limits, as the aperture sizes are greater than our 1000 AU 
aperture. Only the A = 850/im and greater data points resolve IRS 3B from its companion . 



(O'Lineer et al. 


2006) i 


are from 


Barsonv et al. 



and the A =2. 7mm data point is from lLooney et al] (120001 ) . 
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Fig. 8. — SED of of best fit model for L1448-mm A at apertures of 1000 (dashed line) 
and 5000 AU (solid line), plotted with photometry. Data points represented as boxes were 
included in fitting, triangles were not used in fitting but are included for completeness. 
Wavelengths greater than 8.0/mi are upper limits, as the aperture sizes are greater than our 
1000 AU aperture. Only the A = 2.7mm data point resolves the source itself, the rest include 
flux fro m companion sources. A = 12, 25, 60, 100, 350, 450, 800, and 1300/xm data p oints 
are from lBarsony et al.l (119981 ). and the A = 2.7mm data point is from IVolgenaul (12004] ) . 
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Fig. 9a. — Zoomed in images of cavities modeled for each source: (a) IRS 2, (b) IRS 3B, and 
(c) L1448-mm. Images are 2' (30,000 AU) on each side; the primary source being modeled 
is marked with a cross. 
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Fig. 9b. - 
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Fig. 9c. - 
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Fig. 
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10. — Flux versus angle plots of IRS 2A with the best fit model. Models are solid lines, 
are dashed lines. The entire angular range was used for fitting. 
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Fig. 11. — Flux versus angle plots of IRS 3B with the best fit model. Models are solid 
lines, data are dashed lines. The center of the observed cavity is at ~ 1.57 radians. The data 
from to ~ 2.1 radians were used for fitting. Confusion with IRS 3 A and its outflow which 
contaminate ~ 2.1 to 3.14 and to -3.14. 
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Fig. 12. — Flux versus angle plots of L1448-mm A with the best fit model. Models are 
solid lines, data are dashed lines. Only the data from to 3.14 radians were used for fitting 
because data from to -3.14 are contaminated by the binary companion. 
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Fig. 13. — IRAC channel 2 image of L1448-mm (left) compared to the MIPS 24/im image 
(right). L1448-mm A is marked with a cross and L1448-mm B is marked with a star. MIPS 
image courtesy of the cores2disks Spitzer Legacy program. 



